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Abstract. We report on recently performed simulations of Causal Dynamical Triangulations 
(CDT) in 2+1 dimensions aimed at studying its effective dynamics in the continuum limit. Two 
pieces of evidence from completely different measurements are presented suggesting that three- 
dimensional CDT is effectively described by an action with kinetic term given by a modified 
Wheeler-De Witt metric. These observations could strengthen an earlier observed connection 
between CDT and Hoi^ava-Lifshitz gravity. One piece of evidence comes from measurements of 
the modular parameter in CDT simulations with spatial topology of a torus, the other from 
measurements of local metric fluctuations close to a fixed spatial boundary. 



1. Introduction 

Causal dynamical triangulation (CDT) provides an explicit realization of the sum over space- 
time geometries which appears in the path-integral approach to quantum gravity. Its simplicity 
allows us to feed the model into a computer and find expectation values of observables 
numerically. We attempt to exploit the simplicity to address one of the hardest challenges faced 
by non-perturbative approaches to quantum gravity: showing that at low energies a classical 
geometrical theory emerges. 

CDT has provided several pieces of evidence suggesting that it might possess a continuum 
limit showing classical behaviour. However, at present it is unclear what effective theory governs 
this classical regime. In this paper we describe two distinct sets of observables which probe the 
semi-classical regime of CDT in 2+1 dimensions and allow us to obtain non-trivial information 
about the character of the effective theory. Before describing these, let us first briefly introduce 
the set-up of CDT in 2+1 dimensions. 

The motivation for CDT comes from the assumption that path integrals over Lorentzian 
geometries can be computed with the help of statistical path integrals over Euclidean geometries 
that are equipped with a time-foliation. The leaves of a foliated geometry are interpreted as 
spatial geometries and their topology is assumed not to change in time. In CDT the path 
integral is replaced by a discrete sum over all such geometries that can be assembled from a 
finite number of equilateral tetrahedra. More precisely we have the partition function 

ZcBT= e-^cDT[T] ^i^i^ ScMT] = kzN-i-koN^. (1) 

triangulations T 

The action S'cdt is derived from the Regge action and depends linearly on the number 
of tetrahedra and the number A'o of vertices in the triangulation T. The three-dimensional 

^ Based on a talk given at Loops '11 in Madrid, Spain on May 24th, 2011. 



triangulations are required to be foliated in the sense that they can be constructed from a 
sequence of two-dimensional triangulations connected by slices of tetrahedra (see figure |4]) . 

We use Monte Carlo methods to randomly generate CDT configurations according to the 
Boltzmann distribution in ([T]). By performing measurements on individual triangulations we 
can approximate expectation values of observables. A particularly simple observable is the 
volume V{t) of the spatial triangulation at time t. It has been previously observecj^ that its 
expectation value {V{t)) for CDT with spatial topology of a 2-sphere corresponds to that of a 
proper-time slicing of a round 3-sphere, which is the Euclidean counterpart of de Sitter space. 
Both the expectation value and its quantum fluctuations turn out to be well described by an 
effective action 

Scs[V] = Jdt ^co|r - ciV^ with co,ci > 0, (2) 

which is very similar to the (Euclidean) Einstein-Hilbert action J d'^x^{—R + 2h) evaluated on 
a spherical cosmology ds^ = dt^ + V{t)d^'^ . The only difference is an overall minus sign, which 
ensures that for fixed 3-volume ^ is bounded below while the Einstein-Hilbert one is not. 

This is a general property of the Euclidean Einstein-Hilbert action and is known as the 
conformal mode problem. We can illustrate this by writing the metric in proper-time gauge 
ds^ = dt^ + gabdx°'dx^ , in which case the Einstein-Hilbert action takes the form 

S^^ = Kjdtj d^x^g {^^gabQ''''''gcd - + 2a) , g'^'^" = \ {g'^'g'" + g'^'^g'') -g'^'g'^, (3) 

and Q"-^'^'^ is referred to as the Wheeler-De Witt supermetric. This metric on the linear space of 
metric deformations is indefinite: it is positive definite on traceless deformations but negative 
definite on conformal deformations. As a consequence ^ is not a suitable ansatz for an effective 
action for CDT: we want such an effective action to describe, at least in a semi- classical way, 
the quantum fluctuations we observe in our simulations and to do that it needs to have a (local) 
minimum at the classical solutions. 

Inspired by the Horava-Lifshitz approach to gravity [3j and the presence of a preferred time- 
slicing in CDT, we propose a modiflcation of the ansatz (|3]). Namely we modify the Wheeler-De 
Witt metric by inserting a free parameter A, 

gabcd _^ gabcd ^ 1 (^gac gbd ^ gad gbA^ _ Xg^h gCd _ 

This way we obtain the most general ultralocal kinetic term compatible with foliation-preserving 
diffeomorphism invariance. This metric is positive deflnite in the regime A < 1 /2, while in general 
relativity we have the value A = 1. 

In the remainder of this paper we will put this ansatz to test by measuring observables that 
are only sensitive to the kinetic term in the effective action. 



2. CDT on a torus 

We would like to see what consequence our ansatz has for cosmological models. Clearly the 
spatial volume is an observable which only depends on the conformal degrees of freedom in the 
metric. In order to probe the kinetic term in a non-trivial way, we need another observable 
that depends also on traceless degrees of freedom, i.e. an observable that measures the shape 
of space. Taking the spatial topology to be that of the torus, we have a prime candidate for 
such an observable: the complex modular parameter r describing the conformal structure. It is 
well-known that any metric on the torus is conformally related to a unique flat metric of unit 
volume. Up to diffeomorphisms the latter form a 2-parameter family and can be obtained by 
gluing sides of a parallelogram in the Euclidean plane (see figure [T]) . 

^ See [T] for results in 2+1 dimensions and 2! for more recent results in 3+1 dimensions. 




Figure 1. A harmonic embedding of a 
triangulation of the torus. The modular 
parameter r = ri + ir2 is associated with 
the shape of the parallelogram in the manner 
shown. 
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Figure 2. The quantity deduced from 

direct measurement (in black) and from the 
correlations of ri (in orange) and T2 (in red) . 




In [1] an algorithm is described to associate a modular parameter to a triangulation of the 
torus. Basically, this is done by constructing a periodic embedding of the torus in the Euclidean 
plane which we require to be harmonic, by which we mean that any vertex is located at the 
centre of mass of its neighbours. Such an embedding is unique up to scaling, rotations, and 
translations if we insist that the sum of all squared edge lengths is minimized (figure 1). 

Let's see how this modulus appears in the effective action according to our ansatz. Restricting 
the kinetic term given by the modified Wheeler-DeWitt metric Gx to the spatial volume V and 
the modulus r = ri + iT2 we obtain 

In a CDT simulation we can perform direct measurements of the functional A[g] on the spatial 
triangulation at time t, the result of which is shown in black in figure [2] for several different 
boundary conditions (which we will not discuss here). However, we can also deduce the ratio 
(1/2 — A)^[5] of the prefactors in the kinetic term by measuring the correlation functions 
{V{t)V{t + At)) and {Ti[t)Tj{t + At)) for small At. The results are shown in orange and red in 
figure [2| where we have determined (for all four graphs simultaneously) the overall factor 1/2 — A 
by a best fit. We see that the overlap is reasonably good, which is in favour of our ansatz. This 
way we have determined A for various values of the CDT coupling as is shown in blue in 
figure |3) 

3. Extrinsic curvature at a fixed boundary 

Another more local way of probing the kinetic term is by studying the geometry near a fixed 
boundary. It can be shown that the fiuctuations of the spatial metric Qabi^t) a small time At 
after a fixed initial two-dimensional surface are governed by the inverse ^^^^^ of the modified 
Wheeler-De Witt metric in Q, 

(5a6(At,x)5cd(At,y)) - (5„fe(At,x))(5cd(At,y)) cc 5{x - y)gi,a- (6) 



The natural way to test this ansatz would be to consider three-dimensional CDT configurations 
of which we fix the first spatial triangulation and observe how the geometry of the second 



Figure 3. Determined values of A from the 
minisuperspace method (in blue) and from 
extrinsic curvature (in red). 



Figure 4. The extrinsic curvature at an edge 
in the initial boundary is proportional to the 
number of tetrahedra connecting to it. 



spatial triangulation fluctuates. However, defining observables that capture these fluctuations 
and relating them to the continuum metric turns out to be difficult and highly ambiguous. 
Therefore, we make the assumption that the geometry of the second spatial triangulation is to 
large extent captured indirectly by the extrinsic curvature at the initial surface. Accordingly we 
test the ansatz (|6| in which we replace gati^t) by the extrinsic curvature K^b, but have to keep 
in mind that this substitution might have a systematic effect on A. 

Viewing the CDT configuration as a piece-wise linear manifold, the extrinsic curvature has 
support on the edges and is (up to a constant) proportional to the number A^(e) of tetrahedra 
connecting to the edge e (figure |4]). We have checked numerically that the fluctuations in A'^(e) 
only correlate locally and the correlation does not depend on the system size. This means that 
if we send the system size to infinity only an ultralocal correlation will survive. We take our 
boundary to be a regular triangulation of the torus, because this allows us to unambiguously 
extract the ultralocal part of the correlations in N(e). We have determined A by comparing 
this to our ansatz for the correlation of the extrinsic curvature. The results are shown in red in 
figure |3j 

4. Discussion 

At an effective level comparing CDT to the Euclidean Einstein-Hilbert action is problematic 
due to the conformal mode problem. However, if we replace its kinetic term by a non-covariant 
but positive-definite one given by Gx for A < 1/2 we obtain a sensible ansatz and we find 
qualitative agreement with the experiments described. These findings seem to strengthen an 
earlier observed connection between CDT and Hofava-Lifshitz gravity in [5j. 

The discrepancy in figure |3] between the results of both methods indicates that the 
determination of A is quite subtle and probably sensitive to the precise ensemble of triangulations 
used. However, the overall behaviour is clear: as we approach the phase transition in CDT, 
corresponding to a critical coupling ^ 5.6, A increases to 1/2 at which point Qx becomes 
degenerate. 

A more detailed exposition of the presented measurements will appear soon. 
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